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Abstract 

We find constraints on the generation of super-causal-horizon energy per¬ 
turbations from a smooth initial state, under a simple physical scheme. We 
quantify these constraints by placing the upper limit A c ~ 3.0f/n on the wave¬ 
length at which the power spectrum turns over to k 4 behavior. This means 
that sub-liorizon processes can generate significant power on scales further 
outside the horizon than one might naively expect. The existence of this 
limit may have important implications for the interpretation of the small 
scale power spectrum of the Cosmic Microwave Background. 
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I. INTRODUCTION 


The two best candidate theories explaining the origin of structure in the universe are 
“defects” [1,2] and “inflation” [3]. There are many uncertainties in each theory, and the 
overlap between their predictions is considerable. In spite of this, the predictions of each 
theory are constrained, because any model must satisfy the laws of physics. In particular, 
any model is constrained by causality , which plays a fundamentally different role in either 
scenario. 

This difference can be understood in terms of the evolution of two scales, the Hubble 
distance ( H -1 ) and the causal horizon (g?h). The former is the crucial scale for the evolution 
of density perturbations (cf. Equation 82.9 in [4]). The latter is the largest scale on which 
causal interactions can seed density perturbations. We note that the word “horizon” has 
many uses in cosmology: The causal horizon we refer to here is the absolute limit for causal 
interactions. In other contexts, H~ x can also be considered as a horizon, since it is the 
maximum distance over which causal interactions can occur during one expansion time. 

In defect theories, H~ x is of order g?h at all times. Making the usual assumption of a 
smooth initial state, super-horizon-distance and consequently super-Hubble-distance curva¬ 
ture perturbations must be suppressed. During a period of inflation on the other hand, g?h 
grows to some large value, while H~ x remains constant. This means that g?h H -1 for the 
subsequent history of the universe, and therefore that super-Hubble curvature perturbations 
can exist. (Cf. Figure 8.4 in [5]). 

Since the Hubble distance is the crucial scale for the evolution of density perturbations, 
the presence or absence of super-Hubble curvature perturbations may have profound ob¬ 
servable consequences. If so, it will be possible to distinguish between defects and inflation 
in spite of the lack of preferred models. One way to do this has been proposed by Albrecht 
et al. [6]: It is well known that inflation gives rise to oscillations in the small scale power 
spectrum of the Cosmic Microwave Background [3]. They predict that these oscillations will 
not be present for defect theories, provided energy fluctuations are sufficiently suppressed 
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on super-Hubble scales 1 . It is the issue of quantifying this suppression which we address in 
this paper. 

For generality, we phrase our discussion in terms of the horizon distance g?h- It should 
therefore apply to any theory satisfying the laws of physics which starts from a smooth initial 
state. The key quantity is the “energy” t 00 , where is the stress energy pseudotensor used 
in [8] and [9]. is particularly useful, because it describes the flow of energy and momentum 
in the universe in an intuitive way (cf. Section II) and because t 00 is the coefficient of the 
growing mode of matter perturbations. The power spectrum P(k ) is defined via 

<foo(*)fS,(*')} = (2ir) 3 P(k) 6> 3 \k - k') (1) 

where 6^ denotes the Dirac delta in three dimensions. We use the convention 

/(*) = pyw / d3x e "“- x /<*> < 2 > 

The nature of the constraint on P(k) for scales far outside the causal horizon is well known. 
Several authors [9-13] showed that P(k) must fall as k 4 in this limit. However, knowing 
this large scale limiting behavior is not enough. In addition, we need to know where the k 4 
behavior sets in. We can quantify this in terms of a turnover wave-number k c = yT For any 
power spectrum with white noise behavior on small scales and k 4 behavior on large scales 
we can define this turnover by writing the limits as follows 2 . 

P(k) = 


1 This is an active area of research. Numerical simulations described in [7] suggest that in some 
defect models the suppression of curvature perturbations may alter rather than annihilate these 
oscillations. 

2 k c defined here corresponds directly to that used by Albrecht and Stebbins in [14], since their 

“compensation” factor (^ 1+ ^ k has exactly the same limits. 


C ... k dn 1 
c(±y ...k^du ~ 4 
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More general perturbations need not have white noise behavior on small scales, but k c can 
still be defined in a similar way. 

Given a specific defect model, we should be able to work out a value for k c . In the absence 
of specific models, we may still be able to place a lower limit on it. In this paper, we work 
out a lower limit given one simple scenario for causal perturbations, where matter undergoes 
random displacements at each point in space. Although we do not consider a general case, 
we expect that all other methods of perturbing matter will produce power spectra with the 
same or higher values of k c . 

The paper is set out as follows. In Section II we set up a formalism for generating per¬ 
turbations in a fluid, subject to the laws of causality and energy and momentum continuity. 
In Section III we use this formalism to derive a minimum for k c in our simple scenario, and 
consider the relevance of this limit to more general perturbations. In Section IV we draw 
our conclusions. 


II. GENERAL FORMALISM AND K 4 BEHAVIOR 

We begin with a prescription for generating “causal” fluctuations in the energy. We note 
that the stress-energy pseudotensor satisfies an ordinary (non-covariant) conservation law. 
That is 


d, 


V 1 flU 


= 0 


( 4 ) 


The zeroth component is the energy continuity equation. This implies that changes in energy 
from a smooth initial state can only arise as a result of displacing energy elements. Initially, 
consider a homogeneous held with energy r 0 o(*) = 0 V *. We then move n energy elements 
with weight to; from x to the new positions x + A;(*). Doing this at every point in space 
gives the new energy 


T 00 (x) = j d 3 x' ^TO;h^ 3 ^* — x' — A— MS^^x — x'^j 


(5) 
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where M = J2 m i- We imagine perturbations occurring in this way as it is then possible to 
ensure that they also satisfy other physical laws. At the very least, physics requires 

1. Momentum Continuity: We can generate any r 0 o(*) arising via a process satisfying 

momentum continuity by choosing a set of displacement fields A 4 -(®) satisfying 

N 

^m 8 A 8 (*) = 0 (6) 

i = l 

2. Causality: Causality implies constraints on energy perturbations which can be most 
simply stated by saying that no energy element can move further than the causal 
horizon. That is, 

|A 8 (*)| < g?h (7) 


Substituting from Equation 5 into Equation 2 we obtain 

Tm[k] = (d)W / E (w‘ A '<-> - l) (8) 

For completeness, we demonstrate the well known result [9-13] that k 4 behavior follows 
from these assumptions. For modes well outside the horizon we have k.Ai <C 1. Taylor 
expanding the second exponential we find 


1 /» / n 1 n \ 

roo(k) = {0 w /2 / e ~ tk ' X ( ~ l rn t k.A t (x) - - ^ m t (k.A t (x) f + 0(k 3 )) (9) 

\ Z7r ) ' J k i=1 l i=1 J 

Momentum conservation (Equation 6) ensures that first term in this expansion vanishes. 

The leading term is therefore the one proportional to k 2 } implying that the power spectrum 

must go like k 4 . 

Next, we use the formalism developed here to work out a minimum for the turnover in 
the power spectrum under a simple scheme for perturbing matter. 


III. TURNOVER OF POWER SPECTRUM 


In this section, we consider the simplest scheme for perturbing matter which satisfies the 
laws of physics as stated in Section II. In this scheme, two energy elements with the same 
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weight are displaced equally but in opposite directions from each point in space. Further, 
the displacement vectors of energy elements at each point in space are independent. The 
index i then runs from 1 to 2, with Ai = —A 2 = A. Equation 8 reduces to 

To o(k) = ^ 3/2 J d 3 x e~ lk ' x [cos k.A(x) - 1] (10) 

where A is some constant. We want to work out the power spectrum, defined in Equation 
1. We have 

A 2 r 

To 0 (k)To 0 *(k r ) = -—— / d 3 x d 3 x' e ~ l (k-x-k .x ) |- cos fc'Aix) — 1] [cos k'.A(x r ) — 1] (11) 

( 27 t ) 3 J 

To obtain the expectation value of this quantity, we must perform the functional integration 

(Too(fc)Too*(fc')) = J [dA] P [A] T 00 (k) T 00 *(k r ) (12) 

where P [A] is the probability functional of obtaining a displacement held A. In the case 
where the displacement at each point in space is independent, we can separate P [A] into a 
product of probability distributions for the displacement at each point in space. That is 

P[A] = np(A.) (13) 

i 

where we have discretized space so that the index i takes ay through all points and A; is the 
value of the held at the point ay. Using this discretization, and taking dV to be the volume 
of each spatial element, we get 

(TooityToo*^)) = 

-f- [ n [d 3 A t p (A,-)] Y^(dVfe-^- k '^) [cos k.Aj - 1] [cos k'.A k - 1] (14) 

\ Z7r ) J i j,k 

For the j,kth element in the sum, all but the jth and kth integrals can be carried out 
trivially, giving the value 1. Hence we obtain 

(T„ 0 (fc)T w *(C)} = 

x J d 3 A jP (Aj) J d 3 A k p (A*) e -i(k.x j -k'.x k ) [ cos k Aj _ [ cos k < Ak _ (i5) 
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Separating the sum into parts with j ^ k and j = k, and taking the continuum limit, we 
get 


(Too(fc)Too*(fc')) = / d 3 x / d 3 A p (A) e %k ' x [cos k .A - 1] 


(27r) 3 


(27t) 3 


j d 3 xd 3 x'8 {3 \x - x')e - l{k - a 


— fc / .X / ) 


d 3 Ad 3 A'8^ 3 \A — A r )p (A) p (A 1 ) [cos fc.A — 1] [cos k'.A' — 1] 


The first term vanishes, leaving 


(roo(fc)roo* (k')) = A 2 8^ 3 \k — k ') j d 3 A p 2 (A) [cos k.A — 1]" 


The power spectrum is then just 


P(k) = 


d 3 Ap 2 (A) [cosfc.A —1]" 


We can simplify this expression by assuming that the probability distribution is isotropic, 
which implies that P(k) = P(\k\). Choosing k = (0,0, A;) we find 

A 2 r 

P(k) = / dz V(z) [1 — cos kzf (19) 

(27t) j J 

where 

V{A Z ) = j dA x dA y p 2 ( A) 

roo r 

= / dA r dA^p 2 (A) (20) 

J\A Z \ JO 

and A Z} A r and are the components of A in cylindrical polar coordinates. If no energy 

element can move further than some horizon distance g?h, then p(A) must vanish for | A| > 

dn, and V(z) must vanish for z > d h- We then see that the limiting behavior for P(k) is 

' 

7^3 JdzV(z) ...A;>cAh _1 

P(k) = {2v) V ' (21) 

_ ^3 f fdzz 4 V(z) ...k^du- 1 

Comparing the limiting behavior with the form described in Equation 3, we see that the 
turnover of the power spectrum occurs at 
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_ / jdzv(z) y 
c (ifdzz*P(z)J 

Since the integrand in Equation 20 is positive definite, we see that V(z) must be a constant 
or monotonically decreasing function of x. By inspection of Equation 22 we see that k c is 
minimized if we take V(z) to be a constant for 0 < x < g?h- This corresponds to the following 
choice for the probability distribution p( A): 

p 2 ( A) oc h (1) (A r - d u ) (23) 



This is the extreme case where the energy elements always move as far as they are allowed. 
Substituting V(z) into Equation 19 and carrying out the integration, we find the explicit 
solution 


( 4sin£;dH sin2£;dH 
P(k) oc 1-—-h 


(24) 


3kdu 6kdu ) 

By examining the limiting behavior, we find that the value of k c is 204g? h -1 ~ 2.1g?h _1 - The 
corresponding wavelength is A c ~ 3.0g?h- 

In Figure 1 we compare this power spectrum to the results of a computer simulation where 
we perturb matter on a 128 X 128 X 128 lattice by moving pairs of energy elements randomly 
from each point through a distance g?h- We note that there is a very good agreement 
between the theoretical and the simulated power spectra, up to deviations for large k which 
are finite grid effects. This agreement represents a solid check on our calculation. Further, 
it is interesting to compare the continuum calculation described in this paper with the 
discrete case represented by the computer simulation. The discrete case does not involve 
the unphysical assumption that matter at infinitesimally close points in space undergoes 
displacements in arbitrarily different directions. Figure 1 illustrates that this only affects 
the small scale region of the power spectrum, and not its turnover. 

In this section, we have worked out limits on the generation of energy perturbations, given 
minimal assumptions about causality and energy and momentum continuity. Although we 
have only considered special perturbations, we believe that the value of k c we have found is 
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a general minimum. To show why, we consider the two ways in which we could make our 
perturbations more general. 

• Firstly, we could allow more than two energy elements, or energy elements of different 
weight, to be displaced from each point in space. Computer simulations suggest that 
this leaves k c unchanged or slightly increased. 

• Secondly, we could introduce correlations between the displacements of energy elements 
at different points in space. Without performing any detailed calculations, we make 
the following observation. To introduce correlations, we must transport information 
between points which are separated in space. This transport of information takes time, 
which must be subtracted from the time available for energy elements to move. We 
expect that the net effect will be to increase k c . 

IV. CONCLUSIONS 

We have found constraints on super-horizon energy perturbations which can be generated 
from a smooth initial state, under a simple physical scheme. We have quantified these 
constraints in terms of an extremal power spectrum (Equation 24, Figure 1). Its turnover to 
k 4 behavior occurs at a wavelength A c ~ 3.0g?h or equivalently, at a wavenumber k c ~ 2.1 dfi 1 . 
Hence we find that sub-Horizon processes can generate significant power on scales further 
outside the horizon than one might naively expect. Although we have only considered special 
perturbations, we have stated reasons for believing that this limit is general. 

For any theory, constraints on super-Hubble perturbations can be inferred from con¬ 
straints on super-horizon perturbations, provided we know how H~ x and g?h evolve. Es¬ 
sentially, super-Hubble energy perturbations are suppressed for defect theories, and not for 
inflation, but the details will be model-dependent. Since the Hubble distance is the crucial 
scale for the evolution of density perturbations, any difference in this super-Hubble behavior 
may have profound observable consequences, for instance in the small scale power spectrum 
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of the Cosmic microwave background [6]. The limit we have found in this paper makes it 
possible to quantify this without referring to specific models. Therefore, it may be possible 
to discover whether structure in the universe was seeded by defects or inflation, in spite of 
considerable uncertainties in both theories. 
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FIGURES 


Log (P (k) ) 



FIG. 1. Theoretical (solid line) and simulated (data points) extremal power spectrum. Power 
is measured in arbitrary units and k in units of The vertical line corresponds to a mode which 
exactly fills the horizon (kdn = 2ir). 
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